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Abstract 

We study the curvaton scenario using the MSSM flat directions in the gauge-mediated 
SUSY breaking model. We find that the fluctuations in the both radial and phase directions 
can be responsible for the density perturbations in the universe through the curvaton mech- 
anism. Although it has been considered difficult to have a successful curvaton scenario with 
the use of those flat directions, it is overcome by taking account of the finite temperature 
effects, which induce a negative thermal logarithmic term in the effective potential of the flat 
direction. 

1 Introduction 

The curvaton is a scalar field whose fluctuations are isocurvature-like during inflation, but later 
converted to an adiabatic counterpart to explain the structures of the universe ^ El El- The 
amplitude of density fluctuations of the curvaton during inflation is estimated as i3*/<^*, provided 
that the curvature scale is small enough, i.e., V"{4>*) ^ H* > where ff» and 4>* denote the Hubble 
parameter and the amplitude of the curvaton, respectively. The subscript '*' means that the 
variable is evaluated when the cosmological scale is departing the horizon during inflation. 

Among realistic particle theories, the minimal supcrsymmetric standard model (MSSM) is the 
most attractive one, in which there are many scalar fields called fiat directions. The potential of 
flat directions vanishes in the supersymmetry (SUSY) exact limit, but it is lifted by the SUSY 
breaking effect and some nonrenormalizable operators (W ~ (j> n , n — 4 — 9) The flat direction 
(FD) field could be a good candidate of a curvaton. The earlier studies revealed that the curvaton 
scenario is usually difficult to achieve along the context of the gravity-mediated SUSY breaking 
models [21 E]. This is because thermal effects, especially corrections to the potential, lead to the 
fact that the energy density of the FD field will not dominate the universe before it decays. Even if 
thermal effects are suppressed as in the case of the hidden radiation 013], the energy density of the 
FD field does not dominate the universe except for the n — 9 direction, and the amplitude of the 
fluctuations generated during inflation usually damps considerably in the course of the evolution 
after inflation. The crucial point is that the oscillation of the FD field starts rather early so that 
it decays well before it dominates the energy of the universe. 

One may wonder if the oscillation could be 'delayed' by some mechanism. This is exactly 
achieved in the context of the gauge-mediated SUSY breaking models. For many flat directions, 
the sign of the two-loop thermal correction to the potential is negative, which traps the field at 
a large amplitude until very late epoch 0. The energy density of the FD field dominates the 
universe soon after the FD field oscillation starts. The field naturally deforms into Q balls, which 
act as a protector from thermal scatterings, and have relatively long lifetime. In this way, it is 
possible to have successful curvaton models with the use of the MSSM FD in the gauge-mediated 
SUSY breaking models. We will show how it is actually realized below. Notice that the key to the 
problem is the trap due to this two-loop thermal correction to the potential. If it were not there 



(or if it were ineffective) , the energy density of the Q balls would not be able to dominate before 
the big bang nucleosynthesis (BBN) and/or the decay would occur before its domination. 



2 Flat direction 

The scalar potential vanishes along the flat direction $ = 0e ie /\/2, and it is only lifted by the 
gauge-mediated SUSY breaking effect and some nonrenormalizable operators. It can be written as 

® 

where Ms is the messenger scale, Mp — (m^Ms) 1 ' 2 is the SUSY breaking scale, ~ TeV is the 
mass scale of squarks, and Mp = 2.4 x 10 18 GeV is the Planck mass. 

For the curvaton mechanism to work, the mass of the FD field during inflation should be 
negligible compared with the Hubble parameter. It can be achieved in, say, supersymmctric D- 
term inflation models [7] or the no-scale type inflation but here we just assume that there is no 
Hubblc-induced mass term during inflation. Then the FD field slow rolls on the nonrenormalizable 
potential Vnr, obeying the equation of motion 

3£T0 + y£ K (0)~O. (2) 

After inflation, inflaton oscillates around the minimum of its potential, and the universe is dom- 
inated by the inflaton oscillation energy, which behaves as nonrelativistic matter. During this 
stage, inflaton gradually decays into light degrees of freedom, forming dilute radiation. Although 
the energy density of this radiation is very small compared with the total energy density, it actually 
affects the potential of the FD field. There are thermal effects at both one-loop and two-loop order. 
The former one is a thermal mass term which is effective for relatively small amplitude of the flat 
direction : <f> < f~ 1 T 1 where / is gauge or Yukawa coupling constant. The latter one is a two-loop 
thermal logarithmic potential that dominates over the thermal mass term for larger <f>. Hereafter 
we concentrate on the two-loop thermal effects, which can be written as P 1177115] 

V T = c T fTHog(^Pj, (3) 

where ct is a constant of order unity. We expect / ~ 0.1, since the contribution of the gauge 
interactions dominate over that of Yukawa interactions. We found that the sign of ct is negative 
for most of the flat directions, so we assume this is the case, setting ct — — 1. For the sake of 
completeness, we will discuss the positive case in App. A, in which there is no successful curvaton 
scenario. 

Because of this negative thermal logarithmic potential, the FD field will soon be trapped at 
the amplitude 

i 
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It is not released until the zero temperature potential overcomes the thermal counterpart, fT < 
Mp. Notice that the Hubble parameter becomes much smaller than the curvature of the potential 
at that time. 

After the oscillation of the FD field commence, it feels spatial instabilities, and deforms into Q 
balls right after the oscillation begins jSJ 111! 1121 1131 |S] . The oscillation starts at the amplitude 



so that the charge of the produced Q ball becomes 



where (3 < 0.1. * Here and hereafter we assume that tfi osc > Ms ~ Mp/m^ is satisfied. 

If the universe has already become radiation-dominated when T ~ f~ l Mp, the energy density 
of the Q balls the FD condensate is only / 4 times smaller than that of the radiation. Hence the 
energy of the Q balls will dominate the universe soon because it decreases as oc a -3 oc T 3 , while 
the radiation density decreases as oc a -4 oc T 4 . Thus, the energy density of the Q balls becomes 
equal to that of radiation when 

T = T eq ~ fM F . (7) 

On the other hand, if the FD field starts its oscillation during the inflaton oscillation dominated 
(IOD) universe, the energy density of the Q balls evolves as a -3 , while the radiation decreases as 
a — 3/2 rpj^ e energy density of the Q balls will not dominate the universe until 

Notice that Trh f Mp in this case. For the curvaton scenario to work, the temperature T eq 
must be higher than the decay temperature T d (see Eas. (|12|) and l |13f l). 

The Q balls can decay if the mass per unit charge is larger than the mass of the decay particle 
nid- This condition is expressed as 

M F Q~* > m d . (9) 

In the case of the Q ball with B ^ 0, nucleons must be in the decay particles, so that md — ~ 
1 GcV, and the condition becomes 



M F > 



(/^A^m^MJT 3 ) 



„_ 3 \ 2(»-2) 



3 x 10 4 GeV for n = 4, 

9 x 10 5 GeV for n = 5, (10) 
5 x 10 6 GeV for n = 6, 

where we set (3 = 0.1 and A = 1. If the Q ball has no baryon number B = but its constituent 
includes squarks, the decay particles are pions; m N should replaced by ~ 0.1 GeV. If the Q 
ball is non-baryonic, it can decay into neutrinos, and essentially there is no condition like Eq.J^J. 
The Q ball decays by loosing its charge through the surface with the rate ^3] 

so that the decay temperature is given by 

5 , „ 3n-8 

T d ~{T Q M P )> ~__/3-W_£ M P . (12) 
V487T V Mp J 

Imposing Td > MeV so that the BBN is successful, we must have Mp larger than 10 3 , 3 x 10 6 
and 10 8 GeV for n = 4, 5, and 6, respectively. For larger n, Q balls cannot decay before the 
BBN time. Notice that there is another constraint from the gauge-mediated SUSY breaking 

t Since the Hubble parameter is much smaller than the curvature of the potential at the onset of the oscillation, 
we expect < 0.1 instead of /? ~ 10 -4 E|. 
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models; Mp < ^ gm^/ 2 Mp : i- e -> Mp < 10 8 GeV for m 3 / 2 < GeV, where g is a gauge coupling. 
For n = 4, more stringent bound appears as Mp < 10 7 GeV, which comes from the condition 
(t> osc > M s ~ Mp/m^. 

On the other hand, if the Q-ball charge is small enough, such that Q < Q c ~ 10 16 (M f /10 7 GcV)- 4 / n 
|13| . Q-balls decay through thermal effects. In this case, the decay temperature can be written as 

TWO^f. (13) 
Notice that this decay process is realized only for n — 4 with Mp > 3 x 10 4 GeV. 



3 Dynamics of fluctuations 

3.1 Fluctuations in the radial direction 

The second important aspect of the curvaton mechanism is to produce enough amount of the 
fluctuations during inflation, i.e., S = H*/<fi* > 10~ 5 . In general, the position of FD field during 
inflation is constrained as V" < H 2 , leading to (/>* < 4> sr = (HMp~ 3 / A) 1 ^™ -2 ). Hereafter we 
assume that 4>* is of the order of 4> sr for a definite discussion |15j . Thus 

^-a(A)"- 3 , (14) 



p 

leading to 

-J M P , (15) 

H* - \-^5^M P . (16) 

As we will see shortly, however, the amplitude of the fluctuations will damp in the course of the 
evolution after inflation. In general, the equations of motion of the homogeneous and fluctuation 
modes are given by 

4> + 3H<j>+V'((f>) = 0, (17) 
S4> + 3H54>+V" (0)50 = 0, (18) 

where the wave number k — > for the superhorizon mode is assumed for the fluctuation. After 
inflation the homogeneous mode of the FD field still (marginally) slow rolls in the nonrenormalizable 
potential during IOD, and the damping factor is estimated as 

(19) 

where the subscript i denotes the initial values. This slow-roll regime will end when the field is 
trapped by the negative thermal logarithmic potential. Since the curvature at the minimum is 



v"(<p m ) 



KfT) 



2(n-2) 



Mp~ 3 



(20) 



and T ~ (T^HMp) 1 / 4 during inflaton oscillat ion dominated universe, the trap starts when 

(\2 f4,(n-2)rp2(n-2) \ ~ 



4 



The amplitude of the field and the temperature at that moment are 

<f>tr ~ R \ P J , (22) 

( A 5/«-2 T »-i M p)" , (23) 



<-tr 



respectively, so that the damping factor during slow roll in the nonrenormalizable potential is 
estimated as 



ASR) _ 

— 



f 4 T RH y fx 

XMp J \S 



(24) 



This result is applicable if the amplitude of the field during inflation, (/>*, is larger than that of 
the minimum determined by the balance of nonrenormalizable and negative thermal logarithmic 
terms in the potential, namely 

Trh < Tg H = \-^}f- 2 6^^M P . (25) 

In the opposite case, the field fast rolls in the negative thermal logarithmic potential, and quickly 
settles down to the minimum. In this course, the amplitude of 5(f)/(f) does not evolve so much. 
Typically, it only decreases an order of magnitude at most. 

Once the field is trapped at the instantaneous minimum <j) m , it is not released until T < f~ 1 Mp- 
During this stage, the amplitude of fluctuation compared to the homogeneous mode decreases as 




(IOD), 
(RD). 



(26) 



Thus, the additional damping factor during the trap is estimated as 

A3 
(?) 



„ \-2WW=r;f — —T R J" Mp n - 1] M p n(n - 1} , (27) 

where we assumed that the release of the FD field from the trap by the negative thermal loga- 
rithmic potential takes place during radiation-dominated universe, i.e., Trh > / _1 Mjr. Thus the 
amplitude of the fluctuations damps by a factor ^ R \^ during both the slow- roll and the trap 
periods. Since ^f R ^ has dependences on the reheating temperature as 

, T RH ( T RH<Tg H ), ^ 

const. (T RH >Tg H ), 
the total damping factor behaves according to 

er^r4 T t {Trh<TCrh) > (29) 

l (T RH > T£ H ). 

At a glance, higher reheating temperature than T£ H seems to give less damping factor, but it will 
be shown shortly that it does not open parameter space for successful scenario. 
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To this end, let us estimate the damping factor when Trh = T^ H . In this case we can set 
Q, S ~ 1. The amplitude of the fluctuation can be estimated as 



Now let us look at each case of n — 4, 5, and 6 when the Q-ball charge is larger than Q c . For 
n = 4, 6 ~ 10~ 5 can be realized if we take A ~ 1, / ~ 0.03, M F - 3 x 10 4 GeV, and 6 ~ 5 x 10~ 4 , 
for example. In this case, T RH ~ T% H ~ 10 15 GeV. Since T d ~ 30 MeV while T eq ~ GeV, the Q 
balls decay after they dominate the energy of the universe, so that one of the conditions for the 
curvaton mechanism to work is satisfied. However, the gravitinos are overproduced in such high 
reheating temperature, and it is necessary to dilute their density by a factor of (Trj? /T3/2) ~ 10 6 , 
where T 3 / 2 ~ 10 9 GeV for TO3/2 ~ GeV, is the upper limit of the reheating temperature in which 
the cosmological gravitino problem is avoidable |16j . Although there is some entropy production 
when Q balls decay, the dilution factor is only ~ 30, and the abundance of gravitinos is still 
unacceptably large.* 

We can analyze n = 5 and 6 cases in similar manners. For n = 5,<5^10~ 5 can be explained 
by setting A ~ 1, / ~ 0.1, Mp ~ 10 8 GeV, and 6 ~ 2 x 10 -4 . However, the reheating temperature 
is very high ~ 10 16 GeV, which is highly speculative. In addition, the entropy production due to 
the decay of Q balls is not enough to dilute the overproduced gravitinos. For n = 6, the parameter 
set of A ~ 1, / ~ 0.1, Mp ~ 10 8 GeV, and 5 ~ 10 -4 lead to the right amount of the fluctuations, 
i.e., 6 ~ 10 -5 , but still very high reheating temperature such as ~ 10 16 GeV is needed. 

To sum up, n = 4 direction could be a curvaton provided that there was some other entropy 
production. The cases of n — 5 and 6 are hopeless, since too high reheating temperature such as 
~ 10 16 GeV is necessary, where the adiabatic fluctuations of the inflaton should dominate among 
others. As for the larger reheating temperature than T^ H case, it cannot be realized in curvaton 
context for n = 5 and 6, because Trh exceeds 10 16 GeV. For n = 4, Trh can be raised an order of 
magnitude at most, but it only makes 3 times less damping effect, and does not essentially change 
our result here. 

However, the situation changes dramatically in the case of Q < Q c for n = 4 directions, in 
which the Q ball decays through thermal effects. The ideal example is QuQd direction § which 
does not carry any baryonic or leptonic charge, thus we do not have to care about the problematic 
baryonic isocurvature fluctuations. Since there is no damping effect during the slow roll regime in 
the case of n — 4, we only have to estimate the damping factor during the trap due to the negative 
thermal log potential. For A ~ 10~ 2 , / ~ 0.2, Mp = 10 7 GeV, and T RH ~ 10 11 GeV, the charge of 
the Q ball becomes ~ 10 16 , and the decay temperature is as high as ~ 10 3 GeV, while the energy 
density of the Q ball dominates the universe when T — T eq ~ 10 5 GeV. Thus, enough entropy will 
be released through the Q-ball decay (dilution factor ~ 100) to evade the gravitino problem. As for 
the amplitude of the fluctuation, we can obtain £ ~ 10~ 5 for 0* ~ 10 _3 M p and ~ 5 x 10" 5 M p . 
Since the decay temperature is higher than the electroweak scale, the baryogenesis may take place 
in the course of the electroweak phase transition. 

Finally we comment on the case that the FD field is released from the trap during IOD (i.e., 
Trh < f~ 1 Mp). In this case it is impossible to adjust parameters to be (f>* ~ (f>t r , hence the damp- 
ing effects is unavoidable. Also, the later energy dominance by the Q balls is difficult to achieve 
while keeping the amplitude of the fluctuations larger than 10 -5 . Thus, there is no successful 
scenario. 

* Since there is no damping effect during the slow-roll regime in the case of n = 4, it seems possible to have lower 
reheating temperature ~ 10 9 GeV such that the gravitino problem can be evaded. However, the coupling constant 
/ should be unnaturally small ~ 0.003 in order to realize T eq < in this case. 

§The negative thermal logarithmic potential for this direction is shown in App. B. 



—3 




(30) 
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3.2 Fluctuations in the phase direction 

The FD field generally has the A-terms in the potential. For the superpotential Wnr ~ \& n / Mp~ 3 , 
they are written as 

where we assumed the vanishing (or negligible) cosmological constant in the vacuum. As is well 
known, such an A-term appears only if the nonrenormalizable superpotential is <p n -type. If it is 
X4> n ~ 1 -type, where x is anv (super)field other than those consist of the flat direction, no A-term 
results. The x0 n_1 -type superpotential appears for all n = 5 flat directions ^3- (Note that this 
fact does not mean that there is no A-term in the case of n = 5; <p n type superpotential with h > 5 
induces the A-terms, although they are suppressed by ($/JWp)" _Ti . In fact, we have found that 
there is no successful scenario in the case of n — 5.) Therefore, we will consider n = 4 and 6 cases 
here. Since the fluctuation of the potential becomes 

AV(6) ~ Am3 / 2 ^ smn fl S 0\ , (32) 

v ' MjT ]osc 

at the onset of the FD field oscillations, the density perturbation is estimated as 

^~^~A-^ (a^l) f^V^sin^l , (33) 
P V((/) osc ) \ M P J \MpJ losc ' v ; 

where V(4> OS c) ~ Mp and Eq.© are used. 

Now we discuss the evolution of the fluctuations. Since the amplitude of the adiabatic fluc- 
tuations is determined by 9 and 59, not the combination like 59/9, we need follow the evolution 
of 9 and 59 in the potential V(i]) ~ Ani^ /2 <l> n ~ 2 V 2 /M^- 3 , where 7/ = 9 or 59. The mass scale of 
the phase direction, mg, remains smaller than the Hubble parameter until well after the radial 
direction <p is trapped at the instantaneous minimum <f> m , and the amplitudes of 9 and 59 will not 
decrease during the slow roll. The damping takes place only at the very late stage, just before 
<j) is released from the trap. This happens when the slow roll condition on the phase direction is 
violated, me > H, and hence 

T<TsrH (A/^-^l^- 1 ^ 1 ) 5 ^ (IOD), ( ^ 

\f2(n-2) m n-l M n+lj 2 " (RD). 

After the temperature drops down to this value, the amplitude 77 (= 9 or 59) will decrease obeying 
the equation 

77 + 3Hf] + 2^77 + 777^77 = 0, (35) 



where we assumed that 5<f> has died out completely. It is easily derived that 

7, J Hi S« T ^S (I ° D) < (36) 
[ H^^y a. T^tj (RD). 

First we consider the case that the trap ends after the reheating, i.e., Tpjj > T sr . Since the 
damping effect is very mild during the radiation-dominated era compared to inflaton-oscillation- 
domination, one can avoid considerable damping of 9 and 59 in this case. In fact, we obtain a 
damping factor 
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considering both 9 and 59. Thus, the amplitude of the fluctuations can be estimated as 



where is the value of 9 during inflation, and 59 ~ iJ»/</>* is used. It is easily seen that successful 
curvaton scenario is achieved only for n — 6. For example, the parameter set of A ~ 1, / ~ 0.05, 
Mp ~ 10 s GeV, m 3 / 2 ~ 1 GeV, 9* ~ 0.3 and 5 ~ 10~ 2 " realizes the desired density perturbation, 
~ 10~ 5 , where the damping factor £g is ~ 0.3. Notice that the disastrous baryonic isocurvature 
problem is naturally avoided because the baryon asymmetry vanishes due to the nonlinear dynamics 
of the Q-ball formation, in spite of the fact that any n — 6 FD field has non-zero B. This happens 
when the oscillation of the FD field starts after H < me in negative thermal logarithmic potential 

M 

For the parameter set exemplified above, the reheating temperature is constrained as Trh > 
T sr ~7x 10 9 GeV. Such high reheating temperature leads to overproduction of the gravitino by 
a factor ~ {Trh /T 3 / 2 ) ~ 7. However, the Q balls dominate the universe later, and enormous 
amount of entropy will be released. Since the domination by the Q balls begins soon after the 
oscillation of the FD field starts, the dilution factor becomes {H eq / Hd) 1 ^ 2 ~ {pMp/Td) ~ 2x 
10 5 , where T eq ~ 600 GeV and T d ~ 3 MeV (see below Eq.J^), which is enough to dilute 
the overproduced gravitinos. Notice that we just require the reheating temperature higher than 
T sr (RD). Of course, there is an upper limit in order not to have gravitino overproduction. Thus, 
7 x 10 9 GeV < T RH < 10 13 GeV is allowed. U 

Now let us consider the cases that the trap ends during IOD. Then the damping factor becomes 



T 

J- ft-r 



a -1 f -4(fi-l) ( Trh \ 

1 \f-iM F J 



Trh V 4( "- 1} 



m» / - 1 Af? +1 



Tti-1Q 

2(n-l)(3n-2) 



(39) 



The damping is tremendous for n = 4, since £g < 10 -10 even for Mp ~ 3 x 10 4 GeV, / ~ 10~ 2 , 
and TRH/{f~ x Mp) < 1. For n = 6, one could obtain the right amount of the fluctuations. It is 
achieved if we set the parameters as, for example, / ~ 0.1, Mp ~ 10 8 GeV, Trh ~ 5 x 10 s GeV, 
and 5 - 10~ 2 . In this case, T eq - 3 x 10 3 GeV and T d - 1 MeV for /3 ~ 0.1. This situation is 
very similar to the above n = 6 case when the trap ends during RD, since the baryon asymmetry 
is not created because of the nonlinear dynamics of Q-ball formation. However, it is different in 
that there is no gravitino overproduction for relatively low reheating temperature. An example of 
rt = 6 is LLe direction, and the derivation of the thermal potential is given in the App. B. 

Finally we must mention that this mechanism of generating the density perturbation using the 
phase direction is independent of any damping effects on the radial direction, so that a negative 
Hubble-induced mass term can appear during and after inflation. Of course, one must avoid 
Hubble-induced A-terms in order that 9 and 59 should not decrease considerably. They appear 
only if there is three-point interaction in nonrenormalizable Kahler potential, and these terms do 
not exist if the inflaton carries some nontrivial charge. 



4 Conclusion 

We have investigated the possibility of curvaton scenario with the use of the MSSM flat directions in 
the gauge-mediated SUSY breaking models, and shown that the both radial and phase directions 
of the flat directions can act as a curvaton. The later energy domination (before the decay) is 

*6 = H t /<j>* is a little larger than that value determined by the condition that the FD field starts slow rolling, 
but it is possible for the FD field to have an order of magnitude smaller amplitude due to its dynamics )15| . 

II If the messenger scale Mg is less than the reheating temperature, the particles in the messenger sector give 
significant contributions to the gravitino production 1181 . which make the constraint more stringent. However, in 
the present case, Ms ~ Mp/rn^, ~ 10 13 GeV > Trh> and hence the messenger contributions are neglected. 
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usually very difficult to achieve in the realistic particle physics. In addition, thermal effects such 
as corrections to the potential and scatterings leading to early decay (or evaporation) of the field 
make the scenario more difficult to be constructed. In our scenario, both of them are overcome 
by the thermal effects and Q balls. The FD field dominates the universe soon after the release 
from the rather long trap by the negative thermal logarithmic potential. Thermal scatterings do 
not affects so much, since the FD field deforms into Q balls soon after the oscillation starts. The 
Q balls in the gauge-mediated SUSY breaking have a long lifetime, and survive from too early 
evaporation due to thermal scatterings, if the charge is large enough. On the other hand, they 
can decay before BBN (and even electroweak phase transition) for small enough charge. We have 
utilized the intermediate region suitable for the curvaton scenario. 

We have shown that the radial component of the flat direction can play a role of curvaton in 
the case of n = 4. For larger n, S<j)/<p damps considerably in the higher power potential, and we 
have found no successful scenario. In addition, the fluctuations in the phase direction do not damp 
so much and can be responsible for the right amount of the adiabatic perturbation of order ~ 10 -5 
in the case of n = 6, provided that there are no Hubble-induced A-terms. Although such terms 
can appear due to nonrenormalizable three-point interactions ~ J$^$ between inflaton I and FD 
field in the Kahler potential, it can be easily forbidden with use of some symmetry principles. 

Furthermore, there is another problem associated with the baryonic isocurvature perturbation. 
Except for three n = 4 directions, flat directions possess non-zero B (and B—L), which usually leads 
to too much baryonic isocurvature perturbation. In the case of the phase directions as a curvaton, 
the vanishing baryon asymmetry is achieved through non-perturbative dynamics of producing Q 
balls 5 , in which the negative thermal logarithmic potential plays the crucial role. However, 
baryogenesis may be still problematic in this case. In order to avoid large baryonic isocurvature 
fluctuations, the baryogenesis must take place after the Q-ball dominates the Universe. In our 
scenario the Q-ball dominated Universe begins at the temperature T eq ~6x 10 2 or 3 x 10 3 GeV, 
depending on the thermal history (see Sec l3.2l) . Moreover, in order for baryon number to survive 
dilution by the entropy production due to the Q-ball decay, large baryon-to-entropy ratio should be 
generated. Unfortunately, at the moment we do not know such effective baryogenesis mechanism 
which works at low energy scales. 

On the other hand, the n = 4 FD fields such as QuLe and QuQd possess no baryon number, 
so that there is no baryonic isocurvature perturbations from the first place in the case of the 
radial direction as curvaton. Also, the decay temperature can be higher than the electroweak 
scale, therefore the baryogenesis may occur through the electroweak phase transition. We can thus 
conclude that these directions are most promising candidate of a curvaton. 

A Positive thermal logarithmic potential 

Here we comment on the case in which the two-loop thermal correction to the potential has a 
positive sign, and show that there is no successful scenario. We have found that the positive 
potential seems to appear only for some n = 4 directions, so we restrict our discussion only to 
n = 4, although some formulae will be given in general n. The FD field starts its oscillation earlier 
in the positive logarithmic potential, and the charge of the produced Q ball is rather small. It is 
estimated as 




(A.l) 



T, 



(JSC 




(IOD) 
(RD), 



(A.2) 
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(A- 1 / 4 ?!^- 2 )", (IOD), (A3) 
PM P , (RD), 

respectively. The decay process is determined by the evaporation (and/or diffusion) for the Q ball 
with small charge, while the Q ball decays at the rate (|llf> when the charge is large enough to 
survive from the evaporation. Then the decay temperature is 

10^, (Q < Q c ), 

r '~< W, «>«.,, ,A - 4) 

487rQ4 y 

where the survival condition is Q > Q c ~ 10 16 (M F /10 7 GeV)" 4 / 11 [H]. As the energy density of 
the Q ball equals to that of radiation when 

7;„~ ; A-^-^(^)"M F , (IOD), (A5) 
/ 3 M F , (RD), 

it is easily seen that T eq > Td does not hold for n = 4 in any case. Thus, there is no working case 
for the curvaton scenario in the positive thermal logarithmic potential. 



B Thermal logarithmic potential of LLe and QuQd flat di- 
rections 

Here we first derive explicitly the negative thermal logarithmic potential arises taking n = 6 Lyh^e-i 
direction for example. In the end, we also show for n = 4 Q\UiQ2dz direction. Generally, the 
free energy depending the gauge coupling is provided by two loop diagrams as V — Ag^(T)T 4 , 
where A is a numerical factor, and a = 1,2, and 3. Usually, g a {T) can be determined from the 
renormalization group equation 

where t = ln(/i), and b a = 11, 1, and -3, respectively for U(l) Y , SU(2) L , and SU{3) C in MSSM. 
However, in the presence of the flat direction, those fields coupled directly to it acquire the mass 
of order f(j>, where / is gauge or Yukawa coupling. Then, the runnings of gauge couplings (values 
of b a ) change at the scale \i = f<f) accordingly. We will write the new values by b a - Thus, the 
difference of the gauge coupling at the scale fi = T < ftp becomes 

^)^ b ^gl(T)\^ lo S ^y (B.2) 

where g a {T)\^,=o is the coupling without the presence of the flat direction. Since the leading term 
of the potential dependent on <p is written as V = 2g a Ag a (tp)T 4 , the thermal logarithmic potential 
is derived as 

Now we have to evaluate A and b a for the presence of the flat direction. We take an explicit 
example of L\Lie-i for simplicity. Because of the flat direction, only Q, u, d, H u , and SU(3) gauge 
fields are in thermal bath. Since the 53-running will not change, we have to see the changes of 
the runnings of (72 and g\. They are easily derived as 62 = 5 and b\ = 6. The only diagrams that 
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can contribute to the free energy of the type g^T A are those four scalar interactions in the D-term. 
Counting the number of degrees of freedom, we obtain A2 — 5/96 and A\ = 7/288. Thus, 

V T2 U) ~ (^—R--^—]g\ T )THog(^], (B.4) 
- -feffT 4 log (B.5) 

where R — (g' /g) 4 , and the last line is evaluated at [i = T = 10 10 GeV (see Eq.||23)), so / e // w 0.1. 
This is the ideal example for the scenario of the phase direction as a curvaton where the trap is 
released during IOD. 

For Q\U2Qidz direction, we only show the final result: 

M*>*-^(f + !*-^*)rfcn*-*(£), <b,) 

where R 2 — (cte/sa) 4 and R\ = (g'/gs) 4 - Evaluating at \i ~ 2 x 10 11 GeV, we obtain f e ff » 0.2, 
which gives the perfect example for the scenario of the radial direction as a curvaton when Q < Q c . 
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